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Abstract 

We consider the dynamics of a gas of free bosons within a semi-classical Fokker-Planck equation 

c/3 ' for which we give a physical justification. In this context, we find a striking similarity between the 

^ • Bose-Einstein condensation in the canonical ensemble, and the gravitational collapse of a gas of 
g. 

I . classical self-gravitating Brownian particles. The paper is mainly devoted to the complete study 

^ I of the Bose-Einstein "collapse" within this model. We find that at the Bose-Einstein condensation 

^ ■ temperature Tc, the chemical potential ^i{t) vanishes exponentially with a universal rate that we 



CN . compute exactly. Below Tc, we show analytically that \/jMt) vanishes linearly in a finite time 

I tcoii- After tcoii, the mass of the condensate grows linearly with time and saturates exponentially 

' to its equilibrium value for large time. We also give analytical results for the density scaling 

jy-^ ' functions, for the corrections to scaling, and for the exponential relaxation time. Finally, we find 

^ that the equilibration time (above Tc) and the collapse time tcoii (below Tc), both behave like 
— m |r — Tc\, near Tc- 

o 
o 

> 
X 
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I. INTRODUCTION 

The Bose- Einstein condensation is a fundamental result of quantum statistics Below 
a critical temperature Tc, a finite fraction of bosons enters the ground state and a Bose- 
Einstein condensate (BEC) forms 121. Observation of Bose- Einstein condensation in cold 

nnn 

atomic samples has been first reported in 1995 by three different groups [3|, |5| in a vapor 
of spin-polarized ^^Rb, ''Li and ^^Na atoms. Apart from laboratory experiments, another 
interesting application of Bose condensation is related to the problem of boson star formation 



from the dark matter in the Universe Jd| . In that case, Bose star formation involves the axion 
as a dark matter particle candidate [7|. 

Since its discovery 0, 1^, several authors have attempted to develop kinetic models to 
describe the dynamical process of the Bose-Einstein condensation. The dynamical evolu- 
tion of a homogeneous gas of bosons interacting only via "collisions" can be studied using 
an appropriate form of the Boltzmann equation jiol . [ill . Il^ which takes into account the 
specificities of the Bose statistics. Since the gas is isolated and the energy conserved, this 
model describes a microcanonical situation. For this model, the Bose-Einstein condensation 
in .o.entu,„ sp^e Kas been eonsidered by Sen„.o. ^ Tkachev Q and Lac^e et aL Q. 
Alternatively, the canonical evolution of a system of non interacting bosons coupled to a 
thermostat imposing the temperature can be described by a semi-classical Fokker-Planck 
equation Fo. T > r„ wKe.e ,s tKe c.it.ea. te,npe..atn.e, tKis e,nat,on con- 

verges towards the Bose-Einstein distribution. The main purpose of this paper is a complete 
study of this equation below the critical temperature T^, in order to understand the dynam- 
ics of the formation of the condensate. Our study is restricted to a spatially homogeneous 
system without interaction at arbitrary temperature. Alternatively, the collapsing dynamics 
of a trapped Bose-Einstein condensate (BEC) with attractive interaction is often analyzed in 
terms of the Gross-Pitaevskii (GP) or nonlinear Schrodinger equation (see, for instance 
This describes the formation of a spatially inhomogeneous condensate at T = 0. The GP 
equation can display a self-similar collapse, but it occurs in position space while the system 
that we shall consider is spatially homogeneous and the condensation occurs in k-space. 

We thus consider a population of free bosons in d dimensions with dispersion relation 
£:(k) = ^ (we set m = 1 in the following). This system is assumed to be strongly coupled 
to a thermal bath at temperature T. Initially, at t = 0, the system is prepared in an initial 



2 



state which can be, for instance, the equihbrium distribution at some high temperature 
To ^ T. At infinite time, we expect the system to reach thermal equihbrium at temperature 
T, characterized by the Bose-Einstein distribution \IT > Tc (where Tc is the Bose-Einstein 

condensation temperature) 

"^^'"^ = exp(g 

where // is the chemical potential. Below Tc, the equilibrium distribution is the Bose-Einstein 
distribution with = plus a Dirac peak at k = (the condensate) containing the rest of 
the mass. 



A. Bosonic Fokker-Planck Equation 

In this paper, we are interested in the temporal evolution of the occupation number at 
momentum k, that we call p(k, i), from an arbitrary initial state to the final equilibrium state 
described above. If our particles were classical instead of bosons (of course no condensate 
can appear in that case), the coupling to the thermal bath can be modelized by a random 
force and a friction, and the evolution equation for the momentum of a given particle is 
described by the Langevin dynamics 

where i{t) is a delta-correlated random force, whose components satisfy 

{fi{t)f,{t')) = 2D5,,5{t-t'), z,je[l,...,d]. (3) 

In order to recover the equilibrium equipartition theorem, we must impose the Einstein 
relation 

T 

where we have set the Boltzmann constant — 1- The Fokker-Planck equation describing 
the temporal evolution of the momentum distribution reads 



dp_l 



TVi^p + pk 



(5) 



Ultimately, the momentum distribution converges to the expected Boltzmann distribution 

Ps(k) = Z-^expf-^V (6) 



Now, coming back to the Bose gas, we ask whether a stochastic dynamics can be introduced. 



which accurately describes the actual evolution |48|]. Restraining ourselves to non interacting 
bosons, we shall see that one can give a reasonable answer to this question. This will permit 
the description of the dynamical formation of the condensate at and below T^, which is the 
main purpose of this paper. 



Several authors have considered this problem [16|, ^|. Let us introduce a stochastic 
dynamics, by defining the rate at which particles with momentum k get a new momentum 
k'. Following Kaniadakis et al. jl^, we assume the following form of this rate 

Vr(k ^ k') = w(k, k - k') a[p(k, t)] 6[p(k', t)] . (7) 

Contrary to classical dynamics, we assume a dependence of W(k k') with p(k, t) and 
p(k',t). Otherwise, we would simply recover the Fokker-Planck equation Eq. (jSJ, as will be 
shown below. In fact, W(k k') being a rate of departure from the state of momentum 
k, it seems reasonable to assume that it should not depend on the population of this initial 
state. Hence, we set a{p) = 1. Then, one can write the master equation describing the time 
evolution of the particle distribution as 

dp 
dt 

Assuming that the evolution is sufficiently slow, and local, such that the dynamics only 
permits values of k' close to k, one can develop W{h — * k') in powers of k — k' in Eq. (jHJ. 
After some algebra, and proceeding along the line of j|l6|, one obtains a Fokker-Planck- like 
equation 



J [p{k', t)W{k' ^ k) - p(k, t)W{k k')] dk'. (8) 



(9) 



b{p)V^{Dp) + pb{p)3 - Dpb'{p)V^p 

where assuming isotropy, we obtain 

Dik) = ^|^(k,Ak)(Ak)2rf(Ak), (10) 

j(k) = -y"^i;(k, Ak)Akci(Ak). (11) 

We now make the simplification that the diffusion coefficient D{k) is momentum-indepen- 
dent, and that the current J is simply proportional to the particles velocity or momentum 



^(k) = J, (12) 
J(k) = I (13) 
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Going from Eq. ()10|) to Eq. ()12|) is a standard assumption in kinetic theory. It is exactly 
achieved when w(k, Ak) is a symmetric analytic function of (Ak + kAt/,^)/-\/At, where At 
is the discretized time step (ultimately At —>■ 0). This function must also decay fast enough 
so that its two first moments are well defined. In the case of the usual Brownian motion 
based on the Langevin equations this function is simply a Gaussian. Using Eqs. (fT^ 
and (fT^ . Eq. (jH)) becomes 



dp_i 



T(6(p)-p6'(p))VkP + p6(p)k 



(14) 



Note that choosing b{p) = 1, i.e. a transition rate which does not depends on the population 
of the new momentum state, leads to the standard Fokker-Planck equation Eq. (0) for 
classical particles. In the general case, the stationary distribution for k 7^ satisfies 

/I b'ipY 



\p Kp) 

which can be integrated, leading to 



VkP = -/3k, (15) 



Kp) 

where p is an integration constant. In order to recover the Bose- Einstein distribution Eq. (pj, 
we must make the unique choice 

6(p) = l+p. (17) 

We conclude that the classical rate is amplified by the factor h{p) > 1, which translates the 
tendency of bosons to fill already occupied states. Note that the choice 6(p) = 1 — p leads 
to the Fermi-Dirac distribution [l3|, where b{p) now suppresses the probability to hop to 
an already occupied state, as can be expected for fermions (note that a generalized Fokker- 
:'lanck equation incorporating an exclusion principle has been introduced independently in 
in the very different context of the violent relaxation of collisionless stellar systems 
and 2D vortices). This approach has been extended to the case of intermediate statistics, 
associated to b{p) = 1 + Kp (/« G [—1, 1]) 1^. It can be also generalized in order to recover 
Tsallis statistical thermodynamics at equilibrium from a general kinetic equation • 
Focusing on the bosonic case, we finally obtain the associated Fokker-Planck equation 



dp_i 



TVi^p + p(l + p)k 



(18) 
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This bosonic Kramers equation could be directly obtained from a modified Langevin equation 

^ = -^(l+P(k,t)) + f(t), (19) 

where the friction term which tends to move the momentum toward k = now increases 
with the occupation number at momentum k. In this context, the origin of the Bose- Einstein 
instability is quite clear. 

Finally, for later reference, we write the bosonic Fokker-Planck equation for a spherically 
symmetric solution. Setting ,^ = 1, we get: 

For the integrated density 

M{k,t) = f p{k',t)k"^~Uk', (21) 
Jo 

we obtain 

dM (d'^M d-ldM\ ,dM f I DM \ 



dt \ dk"^ k dk J dk \k'^~^ dk 

It has to be noted that the bosonic Fokker-Planck equation (|TH|l valid for non-interacting 

bosons in contact with a heat bath (canonical ensenble) is the counterpart of the bosonic 
Boltzmann equation jl^ 11,Q| valid for an isolated system of interacting bosons (micro- 
canonical ensemble). These kinetic equations are semz-c/asszca/ equations where quantum 
mechanics enters only through the Bose- Einstein statistics. The bosonic Boltzmann equation 
can be derived from a fully quantum treatment as discussed in Sec. IVll The justification 
of the bosonic Fokker-Planck equation from a fully quantum mechanics treatment would be 
interesting but is beyond the scope of this paper. 



B. Analogy with a self-gravitating gas of Brownian particles 

29, 



In a series of recent papers 



281 
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two of the present authors have introduced 
and systematically studied the dynamical properties of a self-gravitating gas of Brownian 
particles in all spatial dimensions. This is the canonical version of the original and certainly 
more challenging problem of self-gravitating Newtonian particles in the microcanonical en- 
semble. In the overdamped limit, the Langevin equation for the position in real space of a 
particle reads 



dr 

'di 



(23) 



where the gravitational potential $(r, t) must be computed self-consistently using the Pois- 
son equation 

A^v,t) = S,Gpiv,t), (24) 

where Sd is the surface of the unit d-dimensional sphere, and G is Newton's constant. The 
associated Fokker-Planck-Poisson (or Smoluchowski-Poisson) system is obtained straightfor- 
wardly 

dp 1, 



at 



TVp + pV$ 



(25) 



A<l> = SdGp. (26) 

As discussed in ji35], these equations also describe the chemotaxis of bacterial populations 
in biology, by a proper re-interpretation of the parameters. From now on, we get rid of non 
essential constants by setting 

G = ^ = Sd = l. (27) 

For a time-dependent radial solution, the system of equations Eqs. ()25|26|1 can be put into 
a unique equation 

^ = t(^ + + p2 + (28) 

dt \ dr"^ r dr I r'^~^ dr ' 



where 



M{r,t)= r p{r')r"^-^dr\ (29) 







is the integrated density. Actually, the equation for M(r, t) looks even simpler 

dM _^ fd'^M d - 1 dM\ M dM 
dt \ dr"^ r dr J r^^^ dr 

It is clear that Eqs. ()28pH()j) are strikingly similar to the dynamical equations found in the 

context of the Bose- Einstein condensation Eqs. ()2()p22|) . Apart from the identical diffusion 

term, we note that the non linear terms have the same dimension as in Eqs. fl20l28|) and 

dimension M'^x[k,r)~'^ in Eqs. ()22I30|) (we shall see later that the -|-1 term at the end of 

Eq. is essentially irrelevant as far as the dynamics of the Bose- Einstein condensation at 

or below is concerned). 

Before summarizing our results concerning the dynamics of a gas of bosons in c? = 3 

described by Eqs. (I2()p22j) . we would like to mention only the main results obtained for a 

self-gravitating Brownian gas obeying Eqs. ()28|30|) . The comparison with this surprisingly 

close model will certainly prove interesting. 



In (i > 3, Eqs. ()28l30p reproduce the known equilibrium profile for T > Tc. For 
T < Tc, a. dynamical instability arises coinciding with the absence of minimum of the 



free energy (thermodynamical instability) 29l l30|. 



In d > 3, for T < Tc, the density develops a scaling profile p(r,t) = pof{r/ro), with 
Po{t) = Tr^'^it) ~ {tcou — t)^^ and /(x) ~ for x +00. Hence, the central 
density diverges in a finite time tcou- Note that / can be calculated analytically in all 
dimensions, and that corrections to scaling have been evaluated 29, 3^. 



In (i > 3, for T < Tc, and for t > tcoii, the central density is swallowed by an emerging 
condensate at r = 0, which grows like Mo(t) ~ (t — tcoiiY^'^^'^ ■ In this post-collapse 
regime, the residual density obeys a backward dynamical scaling, p{r,t) = pog{r/ro), 
with po(^) = Tr'^'^it) ~ (t — tcou)~^ and g{x) ~ for x — > +00. For large time, 
the central condensate saturates exponentially to the total mass, whereas the residual 
density vanishes exponentially with the same rate, which has been calculated by semi- 
classical technics [31]. We found that the relaxation time r above diverges like 
r ~ K^{T — Tc)^^/^, where is known exactly. Approaching from below, the 
collapse time tcou diverges like tcoii ~ K^{Tc — T)^^/^, where again K_ has been 
computed analytically |32j |. 

In (i = 2, at Tc, the density obeys a dynamical scaling with a known scaling function 
/. The central density diverges like po(^) ~ Ci exp(c2-\/t), where ci and C2 are known 



universal constants 3Ci_3J| (in an unbounded domain, the divergence is logarithmic 
instead of exponential 35|). Below Tc, the collapse dynamics occurs after a finite time 
tcou- At tcoiii a fraction T/Tc of the total mass has condensed at r = 0. Using the Virial 
;heorem, we show that all the mass must collapse at r = in the post-collapse regime 
3^. We mention these results in ci = 2, because the dynamics of the Bose- Einstein 
condensation in d = 3 will share some qualitative analogy with this case. 

We thus expect a surprising parallel between the collapse dynamics of a classical self- 
gravitating Brownian gas, which occurs when the kinetic pressure P ~ T is not strong 
enough to counterbalance the gravitational attraction, and a bosonic gas, which forms a 
condensate at k = 0, when the usual Bose-Einstein distribution is not able to accommodate 
for all the mass, even at zero chemical potential (see Section II. A). However, we note one 
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fundamental difference between these two systems. For t +00, strictly below all the 
mass collapses at r = in the gravitational case, whereas the mass of the condensate is 
temperature-dependent in the case of the Bose- Einstein transition. 

C. Summary of results 

In Section II. A and II. B, we briefly review some basic results concerning the static prop- 
erties of the Bose-Einstein condensation, and introduce a simpler model (SM), for which the 
static and dynamical properties can be analytically studied, and which will be equivalent 
in the dense region to the original Bose-Einstein model (BEM) described by Eqs. fl20l22p . 
In section II. C, we show that Eqs. ()20I22|) maximize the rate of dissipation of bosonic free 
energy, and that the dynamical instability exactly coincides with the thermodynamical in- 
stability giving rise to the Bose-Einstein condensation below Tc. 

In Section III, we address the collapse dynamics at and below Tc. In Section III. A, we 
compute the density scaling profile which is found to be independent of the temperature 
T < Tc, up to a temperature-dependent multiplicative factor. In Section III.B, we show 
that at T = Tc the central density (at k = 0) diverges exponentially with time, with a rate 
which can be calculated analytically. We compute the exponentially decaying corrections to 
the final Bose-Einstein distribution with /i = 0. In Section III.C, we show that for T < Tc, 
the chemical potential vanishes in a finite time and that ^2Tfi(t) ~ c(T)(tcoii — t), where 
c(T) is computed analytically close to Tc. We show that tcou diverges logarithmically as 
T approaches Tc from below. In Section III.D, we address the specific collapse dynamics 
obtained strictly at T = 0. 

In Section IV, we study the post-collapse dynamics arising for T < Tc, and for times 
t > tcoii- We show that the condensate mass initially grows linearly, with a slope which 
can be calculated exactly near Tc. For large times, the mass of the condensate saturates 
exponentially fast to its equilibrium value, with a rate analytically known close to T^. 

In Section V, we consider the relaxation time above Tc, which is computed in different 
limits. This relaxation time does not diverge at Tc, but the time after which this relax- 
ation regime occurs (that we call the equilibration time) diverges logarithmically as the 
temperature approaches Tc from above. 
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II. EQUILIBRIUM PROPERTIES AND A SIMPLER MODEL IN A BOX 



A. Bose-Einstein condensation for a gas of free bosons 

Let us briefly repeat the few steps leading to the Bose-Einstein condensation. Considering 
free bosons, the equilibrium occupation number at momentum k is 

PbeO^) = X -• (31) 

exp + /xj - 1 

The total mass is fixed to unity, and the chemical potential /i is defined implicitly by the 
total mass constraint (for simplicity, we set the geometrical factor 3^ — 1, which amounts 
to fixing M ^ Sd=l) 

POO 

M= pBEik)k'^~^ dk = 1, (32) 

^0 

which can be rewritten as the identity 



d-l 



P"'=r S~\ dk. (33) 
•''0 exp + 

Below a certain temperature Tc, this distribution cannot include all the mass. The critical 
temperature corresponds to fj, — 0. Specializing to the case of dimension d — 3 which is the 
focus of this paper (there is no condensation for d < 3), we find 

Pc-[iy\'/'Q, /3c «^ 2.20494... (34) 

where ({x) — J2k=i k~^ is the Riemann ^-function. Below the momentum distribution is 
the Bose-Einstein distribution with n — plus a Dirac peak at k = containing the excess 
of mass Mo, with 



3/2 



37; -T T-Tr 



^'"-'-(j) '-if'- '''' 

Above Tc, we find the following asymptotics which will be useful later: 
and 

3 

H^-lnT, T^+oo. (37) 
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B. A simpler model in a box 



We now introduce a simpler model (hereafter denoted SM), which will share the same 
properties as our original model, as far as the condensation dynamics properties are con- 
cerned. In the kinetic equation Eq. ((201), replace (p + 1) with p, which should be valid in 
the dense region around k = 0, during the formation of the condensate. We obtain 



(38) 



and the following equation for a radial solution 



The integrated density satisfies 



dM „/(92M d-ldM\ 1 fdM^^ 



dt ^[dk^ k dk )^ kd-^ \dk)' 

which is even more similar to the dynamical equation Eq. (jHUI) for a gas of self-gravitating 
Brownian particles than Eq. ()22j). The stationary solution can be easily calculated 

P(k) = (41) 

2 ^ P 

which amounts to replacing the function expx is the Bose-Einstein distribution by the 
function (1 + x). Note that the distribution ()4H) is the Lorentzian. Since this distribution 
p(k) is not integrable up to infinity in d > 2, we introduce a momentum cut-off 



/c < A. (42) 



As before, the total mass constraint reads 



M = i -^^—dk = l. (43) 

2 + 







From now on, we set c? = 3, as well as in the rest of the paper. Then, the integrated density 
can be explicitly calculated 

k 



M{k) =2T[k- ppT arctan (^-^ J J . (44) 
The density profile Eq. (j4H) can accommodate for all the mass up to the temperature 

/3c = 2A. (45) 
11 



Below Tc, a Dirac peak appears at k = 0, containing the excess mass 

T - T 

Mo = (46) 

Just above T^, the chemical potential vanishes in a manner similar to the original Bose- 
Einstein model (BEM), 

whereas at high temperatures, the two models are qualitatively different, as the particles 
momentum cannot spread up to infinity in the SM. Hence, the chemical potential converges 
at high temperatures 

/i^y, T-^+oo. (48) 
C. Some properties of the bosonic Fokker-Planck equation 

Equation ()18|) belongs to a general class of nonlinear Fokker-Planck (NFP) equations 



considered by Kaniadakis [17|, Frank j36| and Chavanis j37|. These nonlinear Fokker-Planck 
equations arise when the coefficients of diffusion, mobility and friction depend explicitly on 
the distribution function. There are different ways to write these NFP equations. One 



(49) 



convenient form for our present purposes is j3 



^ = -v.j. = v.[ifrvp + -Uv* 



where C{p) is a convex function and ^ > can depend on r and t, e.g. on p(r, t). Here, $(r) 
is a fixed external potential but we can also consider the situation where $(r, t) is generated 
by p(r, t) 3^, like in the case of self-gravitating Brownian particles (see Section ITB]) . Taking 
C{p) = plnp — -^(1 + Kp) ln(l + Kp), we recover the bosonic Fokker-Planck equation (fTHj) for 
K = +1, the fermionic Fokker-Planck equation for k = —1 and the classical Fokker-Planck 
equation for n = (for other values of k, this describes intermediate "quon" statistics). For 
the simpler model Eq. (jHHj) . we have C{p) = — Inp describing logotropes js^. Introducing 
a generalized free energy 

F = E-TS = J p^dr + T J C{p) dr, (50) 

one can show that F{t) is monotonically decreasing: F = — J ^C"(p) J^rfr < 0. Therefore, 
if F is bounded from below, the density p(r, t) will converge, for t +oo, to the stationary 
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solution Peqij) such that J* = (in that case F is called a Lyapunov functional). The 
stationary solution of Eq. is determined by 

Pe,(r) = (C")-M-/5$(r)-«], (51) 

where (3 = 1/T and a is an integration constant. This stationary solution is a critical point 
of free energy Eq. (j^Uj) at fixed mass, i.e. it satisfies 6F+aT6M = 0, where aT is a Lagrange 
multiplier. Moreover, it is shown in j^^l that a stationary solution of Eq. is linearly 
dynamically stable if, and only if, it is a minimum of F at fixed mass. In fact, when $(r) 
is an external potential, the critical points of F at fixed mass are necessarily minima since 
6'^F = J C"{p){5pY dr > 0. Therefore, if a critical point of free energy exists it is the 
only minimum of F and, consequently, F is bounded from below: F[p] > F[peq]. In that 

case, the dynamical equation Eq. will relax towards Peqi'r) for t +oo. ^ 

The NFP equation Eq. can be justified in several different ways 37]. It can be 
obtained from the linear thermodynamics of Onsager by relating the current to the gradient 
of a "chemical potential" a{r,t) = — /3$(r) — C'{p) that is uniform at equilibrium (see 
Eq. (HU), writing J* = l/[P^C"{p)]Va. The current can also be expressed as the functional 
derivative of the free energy, writing J* = —l/[^C"{p)]V{6F/6p) . Alternatively, the NFP 
equation can be obtained by writing dtp = —V ■ J and looking for the optimal current 
J=K which maximizes the rate of free energy dissipation F[3] = J(TC"{p)Vp + V$)J(ir 
for a bounded function Eii[J] = |/ C,C"{p)J^ dr preventing arbitrarily large values of the 
current J. This is the variational version of the linear thermodynamics of Onsager. The 
optimal current, that is solution of 6F + 6Ed = 0, is the one appearing in Eq. ()49|) and 
it satisfies F[J*] = —2Ed[J*] < 0. It indeed leads to the most negative value of F (under 
constraints) since 5^(F + Ed) = S'^E^ = |/ C,C"{p){63y dr > 0. These methods emphasize 
the (generalized) thermodynamical structure of the NFP equation. 

III. COLLAPSE DYNAMICS 

A. General scaling solutions for T <Tc 

We now consider the collapse dynamics occurring when the system is suddenly quenched 
to Tc or below from high temperature. It is clear that both considered models should 



13 



behave in a similar manner near the dense region at k = 0. Let us first consider the 
condensation dynamics of the SM. We write 

M{k,t) = 2T I k - ko aictan I ^]] +F{k,ko), (52) 



, 



where ko{t) is uniquely defined by the condition 



2T 



p{k = 0,t) = -j, (53) 
ensuring that 

F{k 0, ko) ~ k^. (54) 

The normalization condition reads 

M{k = A,t) = 1 = 2T (a - /toarctan (^]] +F(A,fco). (55) 



We now look for a scaling solution of the form 

ko 



F{k, ko) = k^f (^] , (56) 



and introduce the scaling variable 

^ - i (^^) 

This leads to 

M'{k, t) = k'pik, t) = + k^-'f l^-j , (58) 

where, from now on, M'{k,t) will denote the derivative with respect to momentum. In 
analogy with the equilibrium density Eq. (jHJ, we observe that fi{t) = /3/co(t)/2 can be seen 
as an effective chemical potential which is expected to vanish as time increases. 
In (i = 3, the dynamical equation for M{k, t) is 



dt V k J k 
Inserting the scaling ansatz, we find 

= -2koT f arctan(x) - / ) + kk^'^ {af - xf) . 
V + IJ 

It is clear that for a = 1, all the terms of this equation scale the same way, if one chooses 
^0 ~ ^0 ^) i-e. /Jo(^) ~ {tcoii — ty^'^- However, we were able to prove analytically that all 
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the solutions of the resulting scaling equation are non physical, leading to a density going 
to a constant for k ^ ko. In addition, when we wrote Eq. we implicitly assumed that 
the second term of the right-hand side (RHS) is a correction to the first one, which implies 
a > 1. Hence, we now consider a > 1, and matching the leading terms of Eq. (j6(Jj) . we 
conclude that 

k = -ck'^-\ (61) 



and that the scaling function introduced above satisfies 

f" + 2 f^\j ' = 2c (arctan(x) - ) , (62) 

which can be integrated once, leading to 

f(x) = — , \(x^ + 6x^ - 3) arctan(x) + 3x - 2x^ - 4xln (x^ + l)] . (63) 

3 [X'^ + 1)"' LV / V /J 

Note that /'(x)/x^ is exactly the scaling function for the density as can be seen from Eq. ()58p. 
We give below some asymptotic results for /(x) which will prove useful later: 

/(x) = — - — a; + 27rclna; + C(l), a; ^ +00, (64) 
o o 

2c . 8c 7 227c q / ii\ „ .„^x 

= —x^ x^H x^ + Ox^M, x->0. (65) 

15 45 1260 ^ ^ ^ ^ 

Let us now analyze the validity of the scaling regime. We ask that the neglected terms in 
Eq. (jUnj) remain small in the range kg <^ k A. Using the large x asymptotic for /(x), we 
find that the non linear term in / appearing in Eq. (j60|) can be neglected, if 

C"'— «'^o'' A;«A;o-^"-'\ (66) 

implying a > 2. Similarly, the term arising from the time derivative of the residual density 
is negligible provided that 

kko'^ («/ - xf) < k^-^ k<^k^ (a ^ 2), (67) 

^ A; < A (a = 2), (68) 

which implies that a > 3 or a = 2. We will find in the next subsections that the admissible 
values a = 3 and a = 2 are in fact respectively associated to the collapse dynamics at Tc 
and strictly below T^,. 

We would like to emphasize that the present model provides an interesting example where 
one obtains a scaling solution which scales differently from what would have been obtained 
from a naive power counting, assuming that all terms in Eq. (jUUI) scale the same way (which 
would lead to a = 1). 
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B. Collapse at T = 



At Tc, the constant appearing in the scahng function is denoted Cc- Expressing the 
conservation of mass, we find 

1 - 2T, (a - A;oarctan (^Aj j = f{A, ko) ~ k^f , (69) 

~ ^AX-' + 0{kl), (70) 

which imphes that 

a = 3, Ce ~ A-3 ~ T^. (71) 

Inserting this value for a in Eq. (jF)T|l . we find that ko{t) decays exponentially at Tc, with 
a rate Cc- To complete the computation of the scaling function, we need to determine the 
constant Cc appearing as a prefactor in Eq. (jUHj) . and which also controls this exponential 
decay. Although the functional form of the scaling functions are identical for the SM and 
BEM, we will see that the constant Cc will differs for both models, as this constant is not 
entirely determined by the dynamics in the region of high density, where the term (p + 1) 
can be safely replaced by p. To make this point clearer, we now look for a global solution for 
the correction term F{k, ko), valid for momentum much greater than ko, and up to = A: 

F{k, ko) ^ ko(j)ik), ko<t:k<A. (72) 

Note that the large x = ^ asymptotic of the scaling function should match the small k 
behavior of (f){k): 

(Aj \ TTC 
Yj ~fc»fco -^^ok'^- (73) 

For the SM, expanding the integrated density up to leading order in ko 

M{k, t) = 2T, (^k - |A;o) + kom + O (A;^) , (74) 
we immediately find that 



0(A) = ttT,. (75) 
On the other hand, substituting Eq. (f7^ in Eq. pUj) we obtain to leading order in ko'. 

koT, + jA = koim - vrT,) = -cMm - ttTJ, (76) 



dt "^V" k 
so that (j){k) satisfies 



0" + l<j>' - /3eCc (0(A) - 0) = 0. (77) 
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This equation admits the solution 



and expressing the condition of Eq. (f75|l . we determine Cc and the full function (j){k) 

c. = 4.^r^, m-^T^[l-'^^^)- (79) 

It is straightforward to check that obeys the small k behavior of Eq. (|73|) . identical to 
the k ^ kg behavior of the scaling function /. 




k 

FIG. 1: We plot the analytical expressions of (p{k) (see Eq. (I79() . dashed line) for the simplified 
model compared to the numerical solution of Eq. (|85() (full line). We have chosen the size of the 
confining box to be A = 1.79805... so that the small k behavior of the two functions coincides. 

Finally, the collapse at Tc for the SM has been fully analyzed, and the central density 
(jHSl) diverges as 

p(0, t) = 2T^kQ^{t) ~ exp (STT^T^t) . (80) 
For infinite time, one converges to the expected equilibrium solution 

M{k, oo) = 2T^k, p{k, oo) = 2T^k-^. (81) 
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t 



FIG. 2: At T = Tc, we plot p^l'^ik = 0,t) ~ kQ^{t) ~ exp(cct) computed numerically by integrating 
the dynamical equation for M{k,t) (see Eq. ()59p). for the SM (A = 1, top full line) and the BEM 
(bottom full line). In both cases, we find an exponential growth, with a rate in perfect agreement 
with our exact results for Cc (the straight dashed lines in these semi-log plots have a slope equal 
to the theoretical values for Cc). For both models, we obtain an excellent fit of p^^'^{k = 0,t) to 
the functional form p^/'^{k = 0,t) = Aexp(cct) — with A ~ 5.7 and B ~ 0.5 (fit not shown, but 
indistinguishable from data starting from t ~ 0.2). 



Now, this approach can be repeated for the BEM, starting with a similar ansatz 

rk k'2 



M{k,t) 







(a 



dk' + kQ(l){k), 



(82) 



exp yyc^2r 

where ii{t) = Pckl/2 acts like a time-dependent chemical potential. For k ^ ko, we find 
that 

dk' = t f dk' - + O {kt) . (83) 

exp [p^j - 1 



exp 



k'^ 



Using the above result, we can insert the expression of M(fc, t) from Eq. (jH^ in the dynamical 
equation Eq. (|^. and find up to the leading order in k^ 



dM 



koTc 



+ 



Pck 



^tanh(^) 



ki^ik) - vrT,) = -cM^ik) - TTTe), (84) 
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or, after simplification, 
with 

0(oo) = ttT^. (86) 

Now, Cc is selected by imposing that 0(A;) is an increasing function (being the integral of 
a density) and that one gets a fast decay of (pioo) — <f){k) ~ exp{—j3k'^/2). Solving this 
eigenvalue problem numerically, we find 

Cc = 1.38452425... = Cn^T^, C = 1.50380614... (87) 

which leads to 

p(0, t) = 2TJ^f{t) ~ exp {2C^'^T^t) . (88) 

In Fig. ^ we have plotted the function (l){k) for the BEM, compared to its analytical coun- 
terpart for the SM. As expected, they significantly differ only for large momenta. In Fig. |2l 
we numerically confirm the exponential growth of the density of particles at k = for both 
considered models, with the growth rate Cc in perfect agreement with the present theoretical 
analysis. Finally, in Fig. El we illustrate the scaling behavior of the density profile. 



C. Collapse for < T < Tc 

We now repeat the above procedure below T^. For the SM model, expressing conservation 
of mass, we obtain 

1 - 2T (^A - fcoarctan (^^j^ = f{X, ko) ~ k^f (^^^ , (89) 

^ ~ ^AX-' + O (ko) , (90) 
Tc 6 

which implies that 

a = 2, c{T)^n{T,-T). (91) 
For this value of a, Eq. (j6T|l can be integrated, leading to 

ko{t) = c{T){t-tcoii). (92) 
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k/k, 








FIG. 3: At T = Tc, we plot {2Tc)~'^k'^p{k,t) as a function of a; = k/ko for the BEM, for times for 
which the central density p{0,t) = 2TckQ^ is equal to 100x2" (n = 1, 14). We have also plotted 
the corresponding scaling function 1/(1 + x^) (dashed line), which is perfectly superimposed with 
the data collapse. A similar plot for the SM would be indistinguishable, for the scale of momentum 
shown. Below Tc, a similar scaling arises for both models. 

We thus find that, below Tc, a singularity should arise in a finite time, for which the central 
density diverges when t goes to tcoii- 

In order to obtain analytical results, we focus on the collapse dynamics just below T^, so 
that <C 1, hence F{k, ko) <^ 1 and c ^ 1. In this regime, non linear terms in F can be 
safely neglected. For momenta k ^ ko, we define a function (l){k) such that 




(93) 



resulting in the following expression for the integrated mass density 




(94) 



Taking A; = A and matching terms of order O (kg) and O (kg), we deduce that 



Tc-T 
Tc 




c 



(95) 
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Inserting the expression of M{k,t) in the dynamical equation Eq. (jHH), we obtain 

T(4>" + j<p')=7rcT-X(^, (96) 



dt V A; 
or 

0" + 1<P' - f3X (0(A) - 0) = 0. (97) 

We recognize the very same equation for (f){k) as in the preceding section, so that the two 
functions are identical up to a multiplicative constant. Hence, for the SM, we obtain 

A = c. = 4v^^T,^ c{T)=A^UT,-n 0(fc) = ^ (l - '"f^^^f ^ ) , (98) 



and a power law divergence of the central density 

p(0, t) = 2Tkf{t) ~ —— ( ) (t^oii - t)-' . (99) 



1 fTr - '2 



87r2r3 V T, 

For the SM, we can evaluate 0(/c), c(T), and A(T) at the next order in 5 = We 
write 

<P{k) = 6Mk) + S^Mk), (100) 

c(T) = 5co + 5'ci, (101) 

A(r) = \o + 6\i, (102) 

where (poik), cq and Aq = Cc can be easily determined from Eq. ()98|1 . Inserting the above 
ansatz in the dynamical equation leads to a complicated linear equation for (j)i{k). We must 
impose the boundary conditions 0i(O) = 0'i(O) = and 0i(A) = 0, which fix the value of Ci 
and Ai. After some cumbersome but straightforward calculations, we obtain 

sinx 



)l X = TT 



A / X 



ci / TT — X \ 27r cos(x) 

I, +1^^+ 3x ' ' 



4tan(x) \J7r t Jstt t / 

vr / r+°° sin(t) , 11 f+°° sm(t) , ^ 

+ - 3/ — —dt / — —dt 

4 V Jsx t 3 Jx t , 

^ 7r(cos(2x) — 1) 



with 



sinfx) 



= !|?(91n(3)-4-9/;-H^*). (104) 
= CiT^, Ci = 13.51919467982... (105) 
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and 

Ai = nCiTl (106) 

Finally, we obtain 

c(T) = 47rT,(T, - T) + C,{T, -Tf + O ((T - T,f) , (107) 
A(r) = A-kT^ + TtCiT^{T^ -T) + {t,{T - T,)^) , (108) 

and note that the next correction to c(T) and A(T) are both positive. 




J I I L_ 



FIG. 4: For the simplified model, we plot 4>{k) = M{k,tcou) — 2Tk obtained by numerically 
integrating the dynamical equation for M(k,t) at T = 0.9Tc (top full line). We have chosen A = 1 
so that Tc = 1/2 and 0(A) = ^k^^ = 1/10 (see Eq. It is in good agreement with our 



analytical expression of Eq. 1)98(1 (top dashed line), which is strictly valid only very close to Tc- We 
also plot the numerical (bottom full line) and theoretical expression of (l){k) (bottom dashed line), 
for the Bose-Einstein model ((/)(A = cxd) = 1 - 27/10^/2 = 0.1462... ^ |^^). 



Again, this approach can be repeated for the more realistic BEM, by writing 

M{k,t)= dk'+[l + -ko)<f){k). 

"'o exp 1(3 — 2~^j — 1 \ / 



(109) 
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At the leading order in ko, the dynamical equation Eq. leads to 



dM 



T 



+ 



(3k 



tanh(^) 



TTCT — Xt 



(110) 



or more simply to 



+ 



(3k 



.tanh (^) 



k 



(f)' -(3X (0(cx)) -0) =0, 



;iii) 



which is again the same equation as the one found in the previous section, which determines 
the function 0(fc) up to a multiplicative constant. In order to compute this constant, we 
have to match the terms of order O (k^) and O (fcg) for k +oo in Eq. which implies 

that 

9.T-T 

(112) 



oo 



-(!)■ 



2 Tr 



(3\' 



(113) 



Finally, we obtain 



3C 



A = c, = C7^2T^^ c(T) = C"7rT,(T, -T), C" = ^ = 2.2557092066. 



(114) 



and 



p{Q,t) = 2Tk,\t) 



C'H'^T^ V T, 



T,.-T 



-2 



{tcoU — t) ^ . 



;ii5) 



Note that at t = tcou, we find that the density and the integrated density are 



2T 6'(k) 

Pik,t,ou) = ^ + ^, M{k,Uou) = 2Tk + (l){k), (116) 



for the simplified model, and 



1 d'ik] k''^ 



for the more realistic Bose- Einstein model. In addition to the equilibrium density, we obtain 
in both models a singular contribution near k = 

(i)\k) TT 



fc2 3A; 
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c(T), ^ 0. (118) 








0,05 0,1 

K 

FIG. 5: For the simplified model (A = 1, bottom line) and the Bose-Einstein model (top line), 
we plot A;o(^) ~^ ~c(T) as a function of kQ{t) as obtained by integrating the dynamical equation 
for M{k, t) at T = 0.9Tc. The theoretical slopes —c{T) evaluated at the first order in Tc — T 
are indicated by an arrow. We find that the numerical value of c(T) are slightly bigger than 
our first order calculation, consistent with the fact that the next correction is analytically found 
to be positive for the SM (see Eq. (^TJ). For small ko{t) {i.e. close to tcou), we find ko{t) ^ 
-c{T) - Acc ko{t), with A ^ 1.35 for the SM, and A ^ 1.65 for the BEM. 

However, we do not observe the appearance of a Dirac peak at k = 0, implying that the 
evolution of the system is not finished yet. Hence, we expect that after tcou, the mass 
included in 0(fc) will be swallowed at = 0, finally giving rise to the condensate. 

Before addressing this issue, we give an estimate of tcoii, near Tc. If one is very close to 
Tc, the density first grows exponentially like at Tc, before crossing over to the behavior of 
Eq. ()115p . when the time is of order tcou- Matching the two regimes, we obtain 

p(fc = 0, t ~ tcoii) ~ exp {2cctcoii) ~ (^^^) t^oii, (119) 
leading to the rough estimate 

t,„,~T-Mn(^^). (120) 
As expected, tcoii diverges as the temperature approaches Tc. This analysis suggests a global 
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form for fco(t), valid for all temperature 

ko = -c{T)h (^^o) , h{Q) = 1, h{x) ~oo X, (121) 

where c(T) ~ Tc{Tc — T) for both models. If we assume the simple form h{x) = 1 + x 
(numerically we find h{x) = 1 + Ax^ for small see Fig. EJ, which is compatible with the 
known asymptotics of f{x) in Eq. (jl2H) . we obtain 

fco = -Ccfco-c(T), (122) 

This equation has a global solution 

ciT) 

ko{t) = (exp [c,{t,oii - t)] - 1) . (123) 

Cc 

For times t close to tcou, this expression leads to the asymptotic scaling of Eq. (jU^ . whereas 
it also reproduces the early exponential decay with rate Cc, which is expected for T close 
to Tc. Imposing that at t = 0, fco should be of order of a typical momentum of the initial 
condition (for instance k^ A for the SM), we recover exactly the estimate of tcoii obtained 
in Eq. (fT2n|l . 

In Fig. I3J we plot the numerical estimate of </)(A;) obtained by numerically integrating 
Eq. at T = 0.9Tc, and compare it to our theoretical results, which are, in principle, only 
valid very close to T^. In Fig. we plot kQ{t) —>■ —c(T) obtained by numerically integrating 
Eq. and find a fair agreement with our theoretical estimates for c(T) , which are strictly 
valid only close to Tc. 



D. Collapse at T = 

For completeness, we now consider the collapse dynamics at T = 0. Since we are inter- 
ested in the density scaling function, we address this question within the simplified model. 
We consider the dynamical equation for the density rather than for the integrated density 



% = ^P^ + ^kpp'. (124) 



Inserting a scaling ansatz of the form 



p{k,t)=po9(^], (125) 



0, 
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in the dynamical equation Eq. 



we obtain 



Pq9 - Po^-xg' = pI (Sfif^ + 2xgg') . 

1^0 



(126) 



We introduce the parameter a defined by 



Po 

Po 



ko 

-a— = apo. 

ko 



(127) 



Then, the equation for the scahng function becomes 



1 

X ' 



g{a - 3g) 

which can be exactly solved, leading to the following implicit equation for g{x) 



(128) 



2a/3-l 



Cx- 



(129) 



where C is a constant depending on the initial conditions, as was already noticed within the 



study of the gravitational collapse at T = 29|, |30|, |3l| . Now, for small x = k/k^, g{x) must 
be an analytic function, with a small x expansion of the form g{x) = g{0) + g"{0)x^/2 + .... 
Matching the small x behavior of the RHS and LHS of Eq. ()129|1 . we obtain 



and 



6 , , 2 
a = -, c/(0) = -, 



gix) - gix)j = Cx-^l\ 



(130) 



;i31) 



We also obtain the exact expression for the central density, and the time-dependent width 
of the dense core k^if) 



p{k = Q,t) = ^po(t) = \{tcoll-t)-\ ko{t) = {t,M-tf'^. 

Finally, at t = tcoiu the small k behavior of the density becomes universal 



(132) 



M{k,t,oii) k^'^\ p{k,t,oii) k- 



-6/7 



(133) 



in contrast with the behavior obtained for < T < Tc, where we found p{k,tcoii) ~ 2T/c ^. 



26 



IV. POST-COLLAPSE DYNAMICS FOR < T < 



As mentioned in Section III.C, the density profile at t = tcou is not yet equal to the 
equilibrium profile which presents a Dirac peak at k = 0. This means that after tcou, the 
density profile should continue to relax to 

T, -T, 



P(k) 



Tr 



for the simplified model, and 

P(k) = 



-(I) 



3/2 



-6{k) + 2Tk' 



(134) 



(135) 



exp (/3f)-l' 

for the Bose-Einstein model. In order to obtain analytical results in this post-collapse stage, 
we again consider the case of a temperature just below T^, such that 



Tr-T 



< 1. 



(136) 



In this regime, and for the SM, the following ansatz is an exact solution of the dynamical 
equation Eq. 



1 



Mo(t)' 



0(A) 



(137) 



M{k,t) = Mo{t) + 2Tk 

where Mo(t) is the time-dependent mass in the condensate, and where 0(A;) has been cal- 
culated in the collapse regime in the previous section. Inserting this ansatz in the kinetic 
equation, we obtain 



dM 



Mn 



\_m\ 



T 1 



Mr 



;i38) 



0(A); 

where we have neglected quadratic terms in 0(/i;), which is justified near T^. For this equation 
to be compatible with Eq. (j97jl . i.e. the defining equation for (pik), we must have 



Mo = c, (0(A) -Mo(t)) 



(139) 



which can be easily solved, leading to the full time dependence of the condensate mass Mo(t) 



Mo(t) = ^ (l - exp [-47r2r3(t - t.^) 



(140) 



For the Bose-Einstein model, we proceed in a similar manner, and find for T close to Tc 



Mo(t) 



1 - 



3 Tr-T 



\ 3/2- 



2 Tr. 



(1 



(l-exp [-CTT^Ti^it- troll)]) 

-Cr^^Tlit-tr 



exp 
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(141) 
(142) 



For both models, the condensate mass initially grows linearly with time 

Mo{t) = nT,c{T){t - t,ou) ~ T^{Tc - T){t - t,ou), t - t,ou < T'^ 



(143) 



whereas it saturates exponentially fast to its equilibrium value, with a rate Cc equal to the 
relaxation rate found at (see Section III.B and the next section), since one can write 

, Mo(t) 



Mo(oo) 

V. RELAXATION TIME FOR T > 



exp [-c,{t - tcoii)] . (144) 



Finally, in this section, we address the problem of determining the relaxation rate to the 
equilibrium solution above Tg. Writing 

^) = /' (J^ \ , + ^M-t/r), (145) 

•^0 exp ( ^ + u ) - 1 



exp + /i 

and neglecting non linear terms in e(fc), we find that 



e{k) = h{pk), (146) 

where h{x) satisfies 

h" + { , „x --\h' + T~'h = 0. (147) 



tanh (f + f ) 



X 



For T ^ Tc (and hence /i ^ 1), we can replace tanh by unity and the solution can be 
exactly written in term of the Kummer confluent hypergeometric function K 



_ . 1 



Imposing a fast decaying solution, we find 

h{x) = x" exp l^-y J , r(T ^ +oo) = ^. (149) 

For the simplified model, the limit of high temperature is not physically interesting since, 
due to the finite box, the particles cannot spread up to large momenta. 

For Q < T — Tc ^ Tc (and hence <^ 1), the eigenequation Eq. ()147|1 coincide with the 
defining equation for the function 0(/c) introduced in Section III.B. Hence, r(T — > Tc) = c~^, 
using the respective values of Cc for the two models considered. 
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It seems paradoxical that the relaxation time does not diverge at Tc, contrary to what 
is expected for a continuous phase transition. However, r is not the equilibration time, 
which is the typical time r^q needed to reach the equilibrium solution. This time can be 
evaluated by considering that close but above T^, the chemical potential initially decreases 
as pk^ity /2 does at and reaches the order of magnitude of its equilibrium value after a 
time Teq. Thus, we write 

koiTeq,Tc) ~ exp(-Ccreg) ~ ko{T) = \f24jiT, (150) 

Using the expressions for /i near Tg given in Eq. ()3fi|l for the BEM and in Eq. ()47|1 for the 
SM, we obtain 

Teq ~ -T-^ln/i ~ T-^ln (^r^) , (151) 

for both models. Note the similarity of this estimate with the expression of tcou near T^, 
given in Eq. p2()|l . The exponential relaxation controlled by r only starts after a time of 
order Teg ^ r. 



VI. COMPARISON WITH OTHER WORKS 



In this section, we give a short review of classical works concerning the dynamics of 
the Bose-Einstein condensation to show how our results compare with these studies. The 
dynamics of the Bose-Einstein condensation has been described by two apparently different 
kinetic theories. The first kinetic theory is based on a quantum version of the Boltzmann 
equation for the one-particle distribution function p(k, t). This is a semi-classical approach 
which introduces corrections for quantum statistics into the ordinary Boltzmann collision 
term Another description is provided by the time-dependent Gross-Pitaevskii 

equation for the condensate wave function 'ip{yi, t) (order parameter for the Bose condensate) 
in a spatially homogeneous medium js^. This equation does not take into account quantum 
fluctuations, or thermal or irreversible effects, but is valid when a large number of particles 
have condensed. 

One of the flrst solution of these kinetic equations is due to Levich & Yakhot who 
consider a gas of bosons without interaction (i.e. neglecting collisions) in contact with a 
thermal bath of fermions. They provide an analytical solution of the corresponding quan- 
tum Boltzmann equation and flnd that the Bose condensate forms in an inflnite time. Since 
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this result is inconsistent with observations, Stoof j41| argues that the quantum Boltzmann 
equation is not vahd to describe the condensate. He considers an isolated gas of bosons 
in interaction (i.e. with interatomic collisions) and argues that the dynamics of the col- 
lapse follows three steps: (i) an incoherent evolution described by the quantum Boltzmann 
equation (ii) a coherent evolution triggering a phase transition (instability) leading to the 
Bose condensate (iii) a thermalization between the condensate and the non condensed atoms 
interpreted as quasi-particles in the sense of Bogoliubov. Stoof focuses on the coherent evo- 
lution. Using a functional formulation of the Keldysh theory, he derives a time dependent 
Landau- Ginzburg equation for the order parameter ('?/'(x, t)) of the phase transition. He also 
shows that the transition temperature for interacting bosons is larger than for an ideal gas 
and that the nucleation of the condensation is short, contrary to the result of |4Q], except for 
temperatures close to the critical temperature Tc. In a more recent paper, Stoof 43] derives 
a Fokker-Planck equation for the probability distribution of the order parameter {ip{:s.,t)) 
which gives a unified description of both the incoherent (kinetic Boltzmann) and coherent 
(Gross-Pitaevskii) stages of Bose-Einstein condensation. A similar program of unification 
of the two theories is carried out by Gardiner & Zoller j^]. Using a projection of the den- 
sity operator p, they derive a quantum kinetic master equation (QKME) for bosonic atoms 
and recover, in particular limits, the quantum Boltzmann equation and the Gross-Pitaevskii 
equation. 

numerically solve the quantum Boltz- 
mann equation and find the formation of a condensate in finite time and the growth of this 
condensate in a post-collapse regime. This is quite different from the results of Levich & 
Yakhot |j4Qi] which are valid when the system is coupled to a bath of fermions. Therefore, 
the original bosonic Boltzmann kinetic equation (without fermion bath approximation) can 
account for a finite time formation of the condensate. More recently, using an analogy with 
optical turbulence, Connaughton & Pomeau obtain a kinetic equation for the spectral 
particle density p(k, t) directly from the Gross-Pitaevskii equation for the wave function 
ipl^Kjt) when the non-linear term is considered as a perturbation (k is the conjugate of x 
in the Fourier analysis). They show that this kinetic equation has the same form as the 
quantum Boltzmann equation. 



Several authors 



m 



141 S . |3l have investigated self-similar solutions of the quantum 



Boltzmann equation in the form p{e,t) = pQ(t) f {e / eo{t)) , where e = j^. In particular. 
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Semikoz & Tkachev j3] find a pre-collapse regime generating, in a finite time tcoii, a dis- 
tribution function of energies scaling as /(e) ~ e^" with a = 1.24. This profile is slightly 
steeper than the Zakharov profile e^^^^ providing an exact static solution of the quantum 
Boltzmann equation corresponding to a constant mass flux J in momentum space toward the 
condensate. The central density increases as po{t) ~ (tcou — ~ (tcou — and 

becomes infinite at t = tcou- The typical core energy eo(t) ~ (tcou — ~ (tcou — t^'^ 

I j 

goes to zero at t = tcou- Semikoz & Tkachev [13] also investigate the post-collapse regime. 
They find that the energy distribution passes from e^" to e^^ after the singularity and that 
the mass of the condensate grows like no{t) = {t — tcou)^^~'^"^^'^^°'~^^ ~ (t — tcou)^'^'^ just after 
collapse. Lacaze et al. [14J also investigate self-similar solutions and obtain similar results 
with an exponent a = 1.234. 

Our approach is physically different since we consider a gas of non-interacting bosons in 
contact with a heat bath establishing a Bose-Einstein distribution (canonical description) 
while the previous authors consider an isolated system of bosons in interaction (microcanon- 
ical description). However, although the dynamics differs in the details (as expected), the 
phenomenology of the collapse is almost the same and the simplified form of our kinetic 
equation (fTHjl allows for a complete analytical solution of the problem (pre and post col- 



lapse) which is not possible for the quantum Boltzmann equation 13|, |lj, |45|, |46|. This is 
clearly an interest of our model. 



VII. CONCLUSION 



In the present paper, we have considered the condensation in k-space of a homogeneous 
gas of non-interacting bosons in a thermal bath fixing the temperature T. In statistical 
mechanics, this corresponds to a canonical description. We have pointed out the striking 
analogy between the dynamical equation describing the Bose-Einstein condensation in the 
canonical ensemble and the one describing the evolution of a self-gravitating gas of Brownian 
particles (or the chemotaxis of bacterial populations). 

In analogy to the case of a classical gas, we have constructed the dynamical Fokker-Planck 
equation for the momentum distribution. For T < Tc, the system experiences a finite time 
singularity, i.e. the density at k = becomes infinite in a finite time tcou ~ T~^\n (^^St'^- 
However, the "singularity contains no mass" and the condensate is formed during the post- 
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collapse regime for t > tcoii- This dynamical scenario leads, for t +00, to a Dirac peak 
in addition to the Bose-Einstein distribution with zero chemical potential, which is entirely 
consistent with the distribution predicted by Einstein at equilibrium Our canonical 
description permits a complete analytical treatment for the different stages of the process, 
contrary to the microcanonical approach starting from the semi-classical Boltzmann equa- 



tion 



13 



14( 1 ■ where one has to mostly rely on numerical simulations. Although qualitatively 



similar (existence of a finite time pre-coUapse regime for t < tcoii, scaling behavior of the 
density, formation of the condensate in the post-collapse regime...), we note that our re- 
sults differ quantitatively from the ones obtained within the microcanonical approach. For 
instance, kQ{t) = c(T)(tcoii — t) in our canonical model, whereas koit) ~ (tcoii — ty in the 
microcanonical ensemble, where, numerically, one finds 7 ~ 1.07 (the difference is even more 
apparent on other quantities) il^. This should not be surprising, as the same phenomenon 
arises in the study of the collapse dynamics of a self-gravitating gas, which is different in 
the canonical and microcanonical ensembles, though qualitatively similar in many respects 

In this paper, we have assumed that the Bose gas is homogeneous. The inhomogeneous 
situation could be studied semiclassically by introducing a mean-field advective term in the 
kinetic equation such that 

where Q{p) is either the bosonic Fokker-Planck operator of Eq. (jl8|) in the presence of 
a thermal bath, or the bosonic Boltzmann operator when elastic collisions are taken into 
account. An even more general model could take into account both contributions, or include 
additional semiclassical contributions arising from the fact that the operators r and k do 
not commute. 

As a final comment, we may note that there exists analogies between the Bose-Einstein 
condensation and the inverse cascade process in two-dimensional turbulence where energy 
accumulates at large scales {k = 0) to form a macrovortex 0|. 
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